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Prompted by the level of accuracy now being achieved in tests of the unitarity of the CKM matrix, 
we consider the possible modification of the Fermi matrix element for the /3-decay of a neutron, 
including possible in-medium and isospin violating corrections. While the nuclear modifications 
lead to very small corrections once the Behrends-Sirlin-AdemoUo-Gatto theorem is respected, the 
effect of the u — d mass difference on the conclusion concerning Vud is no longer insignificant. Indeed, 
we suggest that the correction to the value of |Vud|'^ + |Kjs|'^ + is at the level of 10~*. 



The unitarity of the Cabbibo-Kobayashi-Maskawa (CKM) matrix, U, is a fundamental aspect of the Standard 
Model and any deviation would represent clear evidence of new physics. The most accurate test of unitarity involves 
the (11) element oiU'^U [IHl], namely 



This remarkable result is primarily a triumph of decades of precise studies of super-allowed nuclear /3-decay, which 
at the present time yields a more accurate value of Vud, namely Vud = 0.97425 ± 0.00022, than studies of the free 
neutron. 

In the light of the remarkable level of accuracy that has been obtained, as well as the improvements anticipated in 
the near future, we feel that it is timely to address once more the issue of two corrections which are little discussed. 
The first issue concerns the effect of the mean scalar field in the nucleus. In Ref. 5 it was shown that the difference 
between the radius of a neutron and a proton would grow with density and it was then claimed that this would lead 
to a violation of the Behrends-Sirlin-AdemoUo-Gatto [5| [7] (BSAG) theorem. This result, which was a consequence 
of an incorrect application of the constraint equations in the bag model, is wrong. Indeed, suppose that the u and d 
masses are equal in a nucleon immersed in an isoscalar scalar field. As and md move apart the BSAG theorem 
must apply whatever the scalar field and one must not obtain a correction to the Fermi matrix element linear in 
ruu — rud- Of course, the coefficient of the term in (m„ — m^)^ may change and we estimate this effect here. It is very 
small. 

On the other hand, when we calculate the change in the Fermi matrix element in a manner consistent with BSAG, 
it turns out that at the physical value of m„ — it is no longer totally negligible in comparison with the current 
errors on Vud- The reduction in the overlap of the u and d wave functions turns out to be of order 2 x 10"^. Finally, 
we combine this calculation with an earlier estimate of the effect of pion mass differences carried out using chiral 
perturbation theory, in order to obtain a total correction. When applied to the unitarity test in Eq. ([T]), the total 
correction amounts to about 1.2 x lO^'*. 

As a final note on the possibility of a genuine correction associated with the nucleon decaying in-medium, we point 
out that the relatively small isovector scalar potential (in the QMC model associated with the S meson) expected in 
nuclei with N Z will serve to effectively increase the value of rrid — rUu by an amount Vs. The latter may be as 
large as 4 or 5 MeV in a heavy nucleus and this could produce a decrease in the apparent value of Vud as large as 
2 X 10~^, depending on the neutron-proton asymmetry of the particular nucleus. While for the present this is below 
the level at which it could be detected it is a systematic error which may need to be accounted for in the near future. 



I. INTRODUCTION 




1^ + \Vus\^ + \Vubf = 0.99995 ± 0.00061 . 



(1) 



II. 



EXPLICIT CALCULATION 



In terms of the quark field isodoublet ip 



{''Pu, i^d), the Fermi operator for /? decay is 
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It is related to the generators (T", a = 1, 2, 3) of the isospin transformations 



-I 



T 



which are conserved by strong interactions in the isospin symmetry limit. In the following we assume that the quark 
mass difference 5m = — to„ is the only source of isospin breaking. So we write the exact Hamiltonian of the strong 
interaction as 

H = Ho + v = Ho + ^ J dfTpT^Tp, 

where Hq is the isospin symmetric part: [HqjT"] = 0. Even in the presence of the symmetry breaking term v the 
isospin generators satisfy the current algebra relations 

and combined with the fact that [H, T"] ~ Sm this leads to the BSAG theorem for the Fermi matrix element: 

F+ =< p\J^+\n >= 1 + OiSm^) . 

We wish to estimate how the overlap defect A = 1 — F'^ changes when the neutron is in the nuclear medium and 
the challenge is to do that in a way consistent with BSAG. When we naively use the bag model to compute , as 
originally done in the QMC model, we have the problem that the bag radius depends on the flavor content. As it 
changes during the n — >■ p transition, there is a lack of orthogonality between the eigenmodes of the initial and final 
bag and this leads to an overlap defect of order Sm. Clearly this is an artefact of the static bag model which, by 
definition, cannot handle a radius which changes with time, as is the case in a /3 de(;ay. 

Since this overlap artefact arises from the small mass difference 5m the correct approch is to start with the isospin 
symmetry limit and use the bag model only to compute the deviation from unity. We then find that this deviation is 
explicitly of order (Sm)"^ as it should be. Higher order deviations can also be computed in this framework but we 
stop at order (Jm)^, which is quite sufficient for our purpose. 

We denote as \Ncmr = ±1/2 >, the eigenstates of Hq with isospin 1/2. 

Ho\Na,mr = ±1/2 >= Ma\Na,mr = ±1/2 > 

The nucleon corresponds to a = 0. The exact proton and neutron states are |iV, > such that 

H\N,mr >= M{mr)\N,mr > 
Let P be the projector onto | A^o, = ±1/2 > and Q = 1 — P. We have 

\N,mr >= ^ (l + + • • ■) \^o,mr > 

Z=l-< No\v—^^-—v\No > +0{v^) 
(Mo - HqY 

independent of m^. 

Using T^\No, -1/2 >= |A^o, 1/2 > we have, up to 0{v^) contributions: 

F+ = <N,\\T+\N,-\ > 

The terms linear in v vanish because J^~^ commutes with Hq and hence with Q. From the isospin structure of the 
pertubation we have < Na,—^\v\NQ,—^ >= — < N^, ^\v\No,^ >. So the contribution from the excited states adds 
to that in Z and we end with 



with 
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Since the overlap defect, 1 — F+, is explicitly of order (5m) ^, we can use the (isospin symmetric) bag to finish the 
calculation. If we denote as uja the energies with the same quantum numbers as the nucleon we get 



E 

Q#0 



< a|7o|0 > 



Wo 



< a|7o|0 >= / dr<j)l,-/o4'o , 



with 



where 4>a are the normalized eigenmodes of the cavity of radius R with quark mass m* = fh — g^a and a is the mean 
scalar field in the medium, as calculated for example in the QMC model [HHTU]. In practice we can set m = 0. Note 
that the nuclear vector fields simply shift the energy scale and cannot change F"*". 

If we denote by fla/R the energy (positive or negative) of the mode a, the eigenmode is 



id.fga{r) 



'Att 



with 



/a (r ) - AA„ - Sin [v/f^2 _ (^*^)2 



R 



ga{r) = Na 



R 



{Vta + m*R)r 



sm 



^ni - {m*Ry 



R 



R 



cos 



~ {m*RY 



R 



and M the normalization constant such that 



The energy is determined by the boundary condition f{R) ~ g{R), that is 



(n^ + m*R) sin \^/ni - [m* Rf 



^ni - {m*Ry 



-^^i - {m*RY cos ^Jni - {m*RY 



In the following we assume that m* R is larger than the critical value —3/2 at which f^o = 3/2. This is not a restriction 
in normal nuclei. The scalar matrix element is then 



< a|7o|0 >= 



r"^ ifafo ~ 9c9o)dr . 



We show the resultant overlap defect as a function of density for nuclear parameters [9] in Fig. [T] The u — d mass 
difference has been set to 5 MeV, which is appropriate to the MIT bag [TT]. The overlap defect changes by only a 
few times 10~^ as the density varies from zero to nuclear matter density. 



III. DISCUSSION 



Although the variation of the Fermi matrix element with density is small, we note from Fig. [T] that the intercept at 
zero density is around 1.94 x 10~^, which is only a factor of 10 smaller than the current error quoted for Vud- While 
such a deviation from unity for the hadronic form factor is carefully handled in the s — ?> u transition, it seems to 
receive little attention in the literature in the d ^ u transition. Indeed, we found only the discussion by Kaiser [T2]. in 
the framework of chiral perturbation theory. That work reported a total reduction from the mass difference between 
charged and neutral pions of about 4 x 10~^. The additional contribution from the difference between the ir^nn 
and TT^pp coupling constants, which is of order 0.4% TSl, produces a negligible defect, well below 10~^. The pionic 
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FIG. 1: Deviation below unity of the vector matrix element for the transition neutron to proton, as a function of the density 
of isoscalar nuclear matter. 

correction of Kaiser and the correction from the valence quark overlap calculated here are independent and hence the 
total correction is 6gv = 6x 10~^. Thus one should increase the value of Vud deduced from the analysis of super-allowed 
Fermi /3-decay by this amount. As noted earlier, this amounts to an increase in the value of + iKisP + 

by around 1.2 x 10~^, which is only a factor of 5 below the error quoted in Eq. ([l]). 

While the absence of a sizeable correction to gv in-medium is reassuring, we note that there is a genuine correction 
which has so far been ignored, involving the isovector mean scalar field. This is often ignored because the coupling 
is so much smaller than that for the isoscalar case, but if the mean isovector scalar field felt by the quark is ^V^, 
the isospin breaking term will be proportional to (md — m„ -I- V^)^. This could potentially double the deviation 
associated with the quark mass difi^erence, depending on the neutron excess and structure of the particular nucleus 
under consideration. It will be interesting to explore this effect quantitatively in future as well as to investigate the 
model dependence of isospin breaking term calculated here within the MIT bag model. 
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